Quantum entanglement distribution using a magnetic field sensor 



Marcus Schaffry, 1 Simon C. Benjamin, 1 - 2 and Yuichiro Matsuzaki 1 

1 Department of Materials, University of Oxford, Parks Road, Oxford 0X1 3PH, United Kingdom 
2 Centre for Quantum Technologies, National University of Singapore, 3 Science Drive 2, Singapore 117543 

(Dated: January 27, 2013) 

Sensors based on crystal defects, especially nitrogen vacancy (NV) centres in nanodiamond, can achieve de- 
tection of single magnetic moments. Here we show that this exquisite control can be utilized to entangle remote 
electronic spins for applications in quantum computing; the mobile sensor provides a 'flying' qubit while the act 
of sensing the local field constitutes a two-qubit projective measurement. Thus the tip mediates entanglement 
between an array of well-separated (and thus well controlled) qubits. Our calculations establish that such a 
device would be remarkably robust against realistic issues such as dephasing and multimodal vibrations in the 
sensor tip. We also provide calculations establishing the feasibility of performing a demonstrator experiment 
with a fixed sensor in the immediate future. 



One possible architecture for a quantum computer is based 
on the idea of distributed quantum information processing 
(QIP) [ 1 — 3 1 . Here, rather than controlling all qubits on one 
localized site the set of qubits is distributed to various spa- 
tially separated sites that are entangled with each other. This 
physical distance confers the benefit of better control of the 
individual qubits. Typically for the entanglement operation, 
an optical setup with photons is proposed [4|. In this Let- 
ter, however we show how one can use the dipole-dipole in- 
teraction between electronic spins in conjunction with optical 
detected magnetic resonance (ODMR) to create entanglement 
between different sites. The NV centre defects in diamond are 
very suitable for ODMR and QIP as these possess a long-lived 
spin triplet electronic ground state with the levels |0) and | ± 1) 
that can be easily initialized with a laser, manipulated with 
microwave pulses and read-out optically [5, 6|. This exquisite 
control enables observation of their coupling to adjacent nu- 
clear spins [7 , 8 1 and the measurement of a nearby nuclear 
spin (5). A very promising application of NV centres is their 
capability to detect the strength of very small magnetic fields 
through an induced Zeeman splitting lfT0HT4l . In the follow- 
ing we will show how this high sensitivity to magnetic fields 
can be used to entangle two remote electronic spins. 

Ultimately our proposal is to move a NV centre sensor be- 
tween remote spins to entangle them. We begin by outlining 
a more simple experimental scenario, which could be tested 
in the immediate future. Suppose we are given two electronic 
spin qubits and we can measure the field of these two spins 
by using a crystal defect in a nanodiamond which is placed in 
the middle of the two qubits (see Fig. If one qubit pro- 
duces a field of strength b z at the site of the NV centre then the 
NV centre experiences either —2b z ,Q, or 2b z depending on the 
spin orientation of the qubits: the field is if the two spins are 
in the state | tl) or | |t) and ±2b z if the two spins are in the 
state | ||) or | tt)- Given a sufficiently large external mag- 
netic field, we can prepare the NV centre in the state |+)nv = 
1 / \/2(|0) + 1 1)) . Over time t, this state collects either a phase 
of ±2b z jl-^yt /h = (0±t or depending on the spin state of the 
qubits, where /Inv denotes the magnetic moment of the NV 
centre. Hence this phase allows us to perform a projective 
two qubit measurement. For example, suppose the two qubits 



are initially also prepared in a |+)i/2 = l/v^(| 4-} + | t)) state 
and we let the NV centre precess for the time T = 7t/(Q+, be- 
fore measuring it in the |±)nv = 

l/V2(|0)±|l))-basis. If 
the measurement results in |+)nv> then the two qubits will 
be in the Bell state 1/V2(| it) + I tl)) and similarly, if the 
outcome is |— )nv then the two qubits will be in the Bell state 
1 /V2(\ It) + | tt))- We therefore implemented a determin- 
istic parity projection. As we will show, for times other than 
the ideal T the procedure still achieves a probabilistic parity 
projection. Such a projection is known to be sufficient to im- 
plement the entanglement needed for full QIP |4|. 

(a) static NV centre < ~~ , 
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(b) flying NV centre 



FIG. 1. (a) An NV centre in a nanocrystal is placed in between two 
adjacent electronic spins qubits. The dipole-dipole interaction be- 
tween the spins and a measurement of the NV centre allows to en- 
tangle the two qubits. (b) An NV centre is initially close to one qubit 
and then it is moved near to the second qubit that is far apart from 
the first qubit. Measuring the NV centre at the location close to the 
second qubit again allows entanglement of the two qubits. 

We now analyse this general idea in a more rigorous way 
to show how robust it is with respect to the translations and 
vibrations of the NV centre. First, we look at a configuration 
where the two qubits are located at the origin and at z = 2A 
and the NV centre which is oriented along the z-direction, is 
placed in the middle of the two qubits (see Fig. The 
interaction between the three particles is given by the dipole- 
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dipole coupling: 



"/>.,, = Gjdtf (3(x Lj ■ S,)(i/j • S ; ) - S,- • Sy) 



(1) 



where S, is the spin-operator of the particle i, Xu is a unit vec- 
tor pointing from spin i to spin j, d{ j is the physical distance 
between the spins i and j and Qj = — SMiMy = C is a con- 
stant fl31 . Here jlio is the magnetic constant, jU^/nv = 2\Ib 
are the magnetic moments of the spins, and /Xb denotes the 
Bohr magneton. In an external magnetic field B in z-direction, 
the whole system can be described by the following Hamilto- 
nian: 

# static = H +H DlP with (2) 
Hq = — jU N vBS Zi NV + ^NV^nv — Mi^z.i — ^iBS Zt 2 (3) 
#DIP =#D.NV1 +#D.NV2 +#D,12 , (4) 

here Dnv = 2.87 GHz is the zero-field splitting (ZFS). We 
transform this Hamiltonian to a rotating frame with respect to 
exp(iHot) and neglect fast oscillating terms originating from 
a sufficiently large external field and the ZFS (rotating wave 
approximation). This gives us: 
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static.RWA 
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£>,app,NVl + #D.app,NV2 " 
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(5) 
(6) 



HD.app.ij — 2Cdj j S z .iS z .j 

By using a standard master equation [16| as follows we can 
now evaluate the negative effect caused by dephasing upon 
our basic proposal: 



p'(t) = -i[H staticmA ,p(t)] 
2 



,/=l,2,NV 



-(Sj, z p(t)sl z -^{ P (t),Sy M }) . (7) 



We assume the three particles to be initially in the state 

Iv ( ) = Iv(o)) = |+}iI+)nv|+)2 • 



(8) 



After time t we measure the NV centre in the |±)Nv-basis. 
The qubits are then, depending on the outcome, either in the 
state p + or p_ . This holds for both with the probability 



p± = (l±exp(-f/r 2 , NV )cos(m/2) 2 )/2 , 



(9) 



where a = 2CA . Within the limits of infinite dephasing 
times p± = | \jf±) (\jf±\ are pure where 



W+)=n+\e 
|V-) = l/V2i 
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and where n + is a normalization factor. Our entanglement op- 
eration depends on the measurement outcome. Hence in order 
to quantify the amount of entanglement between the qubits 
that we obtain after measuring the NV centre at time f, we 
calculate the mean entanglement of formation (EF) ifTTIl . i.e. 



which we then plot with respect to the time of the measure- 
ment in Fig. [2] Note that for QIP we could simply discard 
p + retaining p_ which represents a perfect parity projection 
if we know t. In this plot we confirm the expected oscilla- 
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FIG. 2. (Color online) Static NV centre: mean EF with respect to 
the measurement time. For the non-solid curves: T^NV = 2 ms. The 
distance between the NV centre and the qubits is A = 10 nm. 

tory behaviour due to the phase that the NV centre collects. In 
addition one can show that the distortion of the oscillation is 
a result of the qubit-qubit coupling in Eq. Q and due to the 
relatively small magnitude of these terms that they have negli- 
gible effect on the first maximum of the mean EF. Finally we 
conclude that the effect of dephasing for the first maximum of 
the mean EF is also small when the coherence times are suf- 
ficiently large. We can expect large T% times if we implement 
the qubits for example with Sc@C82 or an NV centre. For NV 
centres coherence times of about 2 ms have already been mea- 
sured at room temperature |[T8l[T9ll . Additionally Sc@C82 has 
an unpaired electron spin that mainly exists on the fullerene 
cage ll20l with a coherence time of 200 ms ET1 l22l . Hence 
our entangling operation is capable of achieving high fidelity 
entanglement. 

Next, we address a more practical consideration. Namely 
that the nanocrystal in which the NV centre is embedded can 
be placed on an AFM-tip which then enables us to move the 
defect very precisely between various sites in a distributed 
quantum computer architecture. This ability of spatial control 
comes with the price of limited controllable vibrational modes 
of the AFM-tip, with typical frequencies varying between 
50 kHz and 500 kHz [23|. To characterize the consequences 
of these vibrations onto our entangling operation we consider 
first the effect of the NV centre oscillating between the two 
qubits in a single mode, i.e. c/nv,i/2(/) = A± 5cos(of + 0); 
where 5 denotes the amplitude, (0 the frequency, and (j) the 
phase of the oscillation. As above, we derive a rotating 
wave approximation Hamiltonian and end up with the time- 
dependent Hamiltonian 

#vib,RWA(0 = -Hb,app,NVl (* ) +#D,app,NV2(*) +Hd,\2 ■ (H) 

We can consider this vibrational scenario as a perturbed static 
case and hence we expect a maximum of the mean EF at 
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around tt / joc| (similar to Fig. |5J. Note that for larger times 
the qubit-qubit interaction and the dephasing become relevant. 
For this reason we define the maximal achievable mean EF M 
to be 

M = M(8,co^) = max (EF(r)) . (12) 

0<f<27r/|ce 

We are interested in how the introduction of many modes with 
various parameters (5,(0,0) affect this maximal achievable 
mean EF, and in particular, for which parameter regime M is 
close to 1. To this purpose, we analyse the effect of many 
modes that share the common parameter p on M and define 



9(p)Mt,p))Mt,p)\dp , (13) 



as the mean density matrix with respect to the distribution 
&{p) of the parameter p; where \\j/(t,p)) is the solution of 
the Schrodinger equation for // v ib,RWA(/ )■ We find an analytic 
approximation of this solution as follows. First we neglect the 
qubit-qubit coupling in //yib.RWA(0 

#vib,RWA,app( f ) = #D,app,NVl (? ) +^D.app,NV2(0 • (14) 

We then transform this Hamiltonian to an interaction pic- 
ture and expand it in a power series with respect to 5(f) = 
5 cos (cot + 0). This gives: 

Hi(t) = e' flstetic • RWAr Hvib.RWA,app^ , ^*' RWA, (15) 
= .1 I ((-iy3^ + 6^+O(«(0 3 ))w«J- 
Finally, we apply time-dependent perturbation theory: 



Pi(t) = Pi(0)-i [///(0,P/(0)]*' 
Jo 



+ H) 2 / / [ff / (0,[#/(O.p/(0)K'<fc / , (16) 

Jo Jo 

where P/(0) = and get an analytic approximation of 

the solution of the Schrodinger equation for Eq. ( fTl) . The first 
order approximation in 5 reads: 

,s , n \ - 6C c sin ((Of + 0) -sin(0) 
pi(t) = f>i(0) + i^8—± ^ ^ 

x[S z<m S z ,i-S zm S Z:2 ,pi(0)] ■ (17) 

This approximation allows us to understand the effect of many 
modes. First, we analyse the effect of many vibrational modes 
with different frequencies. In Fig. [3^ we compare our approx- 
imation with the exact solution for ( fTT| by plotting the max- 
imal achievable mean EF for ~p m (in both cases) with respect 
to the vibration frequency CO for different amplitudes and for a 
particular phase. For the average p m in Eq. ( fT3| l we use a trun- 
cated Normal distribution on [0,°°) with mean CO and standard 
deviation a = 0.01(0. In the static case we have seen that the 
mean EF is maximal when we wait about time %j \ a | before 



we measure the NV centre. First assume the NV centre is ini- 
tially in the middle of the two qubits (<j> — n/2). If the NV 
centre oscillates for 1/2 + k periods (k periods), where k is 
a non-negative integer before it is measured, then the asym- 
metry in the system is maximal (minimal) and we get a low 
(high) mean EF. Obviously for a large k, this effect is less 
pronounced. The minus sign in the commutator in Eq. ( fT7| ) 
describes exactly this asymmetry which decreases for large 
frequencies and thus the perturbative solution ( fT7] i explains 
the characteristics of Fig. [3^. We note that analogous plots 
where we increase and decrease the width of the frequency 
distribution are not very sensitive to the width of the distribu- 
tion which is in agreement with the rate in Eq. ( fTTj i. 
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FIG. 3. (Color online) (a) Maximal achievable mean EF for vibration 
modes where the vibration frequencies (amplitudes) are distributed 
according to a truncated normal distribution on [0, °°); with mean CO 
(8) and standard deviation 0.01(0 (0.015), i.e. p® (p ) with respect 
to the mean (common) vibration frequency a>/2n for different com- 
mon (mean) vibration amplitudes 8 and <j) = j . (b) Maximal achiev- 
able mean EF for vibration modes where the phase is uniformly dis- 
tributed, i.e. with respect to their common vibration amplitude 8 
for different common vibration frequencies. In both plots the solid 
lines are calculated from the numerical solution of the von Neumann 
equation for i? v ib,RWA(f) an d the dashed lines are calculated by using 
the approximations discussed in the text. Here A = 10 nm. 

Second we can repeat this analysis by averaging over the 
amplitude instead of the frequency, i.e. by considering p s . 
Again we see from Eq. ( fT3| that averaging over a reasonably 
narrow truncated Normal distribution on [0,°°) has little ef- 



4 



feet, and hence, we get the same plot as in the analysis before. 

Third we analyse the effect of many vibrational modes with 
different phases. Thus we assume a uniform distribution of 
phases on [0,2?r]. With Eq. ( fT6] l we obtain in second order in 
8 the following density matrix 

pf(t) = p 7 (0) - i6~8 2 t[S z , m S zA + S z , NW S z . 2 ,p,(0)} 
_ 18 C% 2 l-cos(a») 

x &,NvS z ,i — S z .nvS z ,2, [S Zl NvS Z: \ — S z ^yS Z: 2,Pl(0)}] . (18) 

Again we compare this approximation with the exact solution 
for ( fTT| by plotting in Fig. the maximal achievable mean 
EF with respect to the vibration amplitude for different fre- 
quencies and find a good agreement between the two solutions 
especially for small amplitudes. We see that the last term in 
( fT8| l can be seen as a (time-dependent) decoherence rate and 
explains why higher frequencies and smaller amplitudes are 
less disturbing to our entangling operation. 

In summary we conclude that if oscillations of the NV cen- 
tre are not avoidable then we can mitigate the adverse effect 
for our entangling operation by ensuring that the oscillations 
are very fast and of low amplitude. 

Having established that a static NV centre can entangle two 
external spins, and having further determined that the realistic 
issues of dephasing and oscillation do not present fundamen- 
tal difficulties, we are now in a position to consider entangling 
two remote spins via a moving, or 'flying', NV centre. There- 
fore we consider a distance D^> A between the two qubits and 
propose that the NV centre should now fly from one qubit to 
the other, i.e. from z = Zo to z = D — zq with velocity v. Hence 
the NV centre interacts first with qubit 1 and then with qubit 
2 before it is measured at tM = (D — 2zo)/ V (see Fig. [TJ)). 

As above we derive with the time-dependent distances 



^NV,i(0 =tv + zo and<i NV ,2(0 = D — tv — Zo , (19) 
the following rotating frame approximation Hamiltonian 

^flying,RWA(0 = #D,app,NVl (0 +^D.app,NV2(0 , (20) 

where we neglect the qubit-qubit coupling due to the large 
distance between the qubits. 

Again we initialize the three spins in the state |*P,-) = 
(see Eq. |8]l). The Schrodinger equation determines the state 
I^F^m)) when the NV centre travelled from one qubit to the 
other. The measurement of the NV centre in the |±)nv -basis 
projects the qubits either to the state 



= | U> + 1 tt»/V2 or 



(21) 



2e'l 



= N + [e'P 1 1|) + | ft) + y— p (I It) + I U» ) (22) 
each with the probability: 
p fly ,_ = sin 2 (j3/2)/2 and Pfly , + = (3 + cos(j3))/4 , (23) 



where j3 = C ((D — zo)~ 2 — Zq 2 ) /v and is a normalization 
factor. The EF of the first state is 1 and that of the second state 
is given by the binary entropy function H 

E F(p + )^(i + i^)wi.^l^.(24, 

Hence whenever j3 = it mod 2n the entangling operation is 
maximally entangling. Roughly, shorter entangling opera- 
tions achieve more entanglement as they are less affected by 
dephasing. Hence we define the optimal velocity v op t for 
our scheme to be the fastest velocity for which the mean EF 
reaches a maximum i.e. /3 (v op t) = 7t. In Fig. [4] we plot the EFs 
and probabilities around the optimal velocity v opt . 
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FIG. 4. (Color online) (a) Flying NV centre: EFs, mean EF and 
success probabilities as a function of the inverse velocity v~' when 
the NV centre has moved from one qubit to the other. Parameters: 
D = 100 nm,zo = 5 nm. 

In summary, for this Letter we present an entanglement op- 
eration that is based on sensing small magnetic fields. Our 
calculations establish that this operation is fast enough to ne- 
glect dephasing for robust qubits. Moreover, imperfections 
in the scheme caused by vibrations of the NV centre are tiny 
for high frequencies and low amplitudes. Eventually remote 
qubits can be entangled by physically moving the NV cen- 
tre mounted, for example, on an AFM-tip between the two 
qubits - making the entangling operation of high value for a 
distributed quantum computer architecture. 
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